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A  CONTINUOUS  BIN  MODEL  OF  AEROSOL 
BEHAVIOR  IN  THE  MARINE  BOUNDARY  LAYER 

Andreas  K.  Goroch 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory 
Atmospheric  Science  Directorate 
Monterey,  CA  93943-5006 


The  nature  of  aerosol  panicles  at  the  ocean  surface  is  important  to  a  large  number  of  direct  and 
indirect  applications.  Visibility  in  the  optical  wavelengths  is  largely  controlled  by  the  distribution 
of  haze  particles.  The  growth  of  ice  and  rime  on  ship  structures  depends  on  the  concentration 
of  aerosol  particles,  as  well  as  temperature  and  water  vapor  content  of  the  atmosphere.  Less 
directly,  aerosols  play  an  important  role  in  the  transfer  of  mass  and  latent  heat  from  the  ocean 
to  the  atmosphere.  Particularly  at  high  wind  speeds,  the  aerosol  population  can  be  responsible 
for  a  significant  part  of  the  total  surface  flux  of  heat  and  moisture. 

Operational  Navy  applications  require  a  particularly  accurate  description  of  the  panicle  size 
distribution.  Since  the  efficiency  of  scattering  by  particles  strongly  depends  on  the  ratio  of 
incident  wavelength  and  scanerer  dimension,  visible,  infrared,  and  microwave  systems 
performance  are  affected  by  different  pans  of  the  aerosol  size  distribution.  In  particular,  micron 
droplets  attenuate  or  modify’  visible  imagery,  while  larger,  submillimeter  particles  affect 
microwave  wavelengths.  The  microphysics  of  the  aerosol  population  is  therefore  important  to 
the  development  of  operational  applications. 

The  behavior  of  marine  aerosols  has  been  studied  by  several  investigators.  Ling  and  Kao  (1976), 
developed  a  monomodai  size  distribution  model  which  described  aerosol  thermodynamics  and 
dynamics  in  the  several  meters  above  the  ocean.  The  assumption  has  been  that  the  state  v  fables 
(temperature  and  humidity)  are  described  by  diffusion  equations  arising  from  the  temperature  and 
humidity  gradients  respectively.  This  work  has  been  expanded  to  other  more  realistic 
distributions  by  Stramska  (1986),  who  has  used  a  discrete  bin  model  of  the  aerosol  size 
distribution.  Burk  (1984)  has  invoked  higher  order  closure  to  simulate  the  effect  of  aerosol 
behavior  on  the  state  variables  This  model  considered  the  conservation  of  marine  atmosphere 
salt  particles,  which  would  accommodate  to  the  ambient  moisture  and  temperature.  This  model 
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provided  the  clear  advantage  that  the  aerosol  growth  equations  correctly  include  the  solute  of  a 
solution  on  particle  growth  rates.  The  detailed  dynamics  of  individual  particles  has  been  treated 
with  a  Lagrangian  model  of  aerosol  transport  proposed  by  Fairail  and  Edson  (1989).  The  model 
is  currently  in  use  to  investigate  the  anisotropy  of  aerosol  source  function  around  the  wave  crest. 

A  common  difficulty  of  these  models  is  the  correct  description  of  the  detailed  vapor  condensation 
and  growth  on  the  marine  aerosol  particles.  The  microphysics  of  a  droplet  population  have  been 
successfully  described  in  the  context  of  precipitating  clouds  by  Young  (1974).  In  this  model,  the 
aerosol  distribution  is  considered  as  a  piece-wise  linear  function,  modified  by  the  ambient  state 
variables.  This  model  has  the  advantage  of  providing  a  detailed  description  of  the  aerosol  size 
distribution,  as  well  as  being  somewhat  more  realistic  in  the  treatment  of  condensation, 
evaporation,  and  coagulation  in  precipitating  clouds.  The  following  paper  applies  the  piece-wise 
continuous  aerosol  model  to  the  behavior  of  marine  aerosols  near  the  ocean  surface.  The  mar  r 
atmosphere  is  considered  as  three  connected  systems:  state  variables,  continue -s  aero;  >1 
processes,  and  discrete  aerosol  processes.  The  state  variables  of  temperature  and  moisture  behav  e 
according  to  the  accepted  first  order  closure  (K  theory)  assumptions.  The  continuous  aerosol 
processes  consist  of  evaporation  and  condensation.  The  discrete  processes  involve  the  changes 
of  numbers  panicles  either  by  motion  to  an  adjacent  layer,  or  by  injection  of  aerosols  from  wave 
breaking  or  other  mechanisms  at  the  surface.  The  model  provides  a  description  of  the  venical 
distribution  of  aerosols  in  the  turbulent  boundary  layer  within  100  m  of  the  ocean. 

Governing  Model  of  the  Atmospheric  Surface  Layer 

Our  model  describes  the  near  surface  marine  atmosphere  by  a  venical  logarithmic  grid.  At  each 
grid  point  we  maintain  the  values  of  temperature  humidity,  and  aerosol  parameters.  Aerosol 
parameters  considered  by  the  model  are  the  number  and  mass  of  aerosol  panicles  in  groups 
defined  by  a  smallest  and  largest  radius.  The  aerosol  parameters  are  modified  by  a  continuous 
change  model  which  describes  the  vapor  condensation  onto  and  the  evaporation  from  ihe  panicle 
surface,  and  the  discrete  change  model  which  considers  panicle  venical  motion  and  injection 
from  the  surface. 


State  Variables  in  the  Atmospheric  Surface  Layer 


The  temperature  and  humidity  behavior  of  the  atmosphere  within  100  m  of  the  ocean  surface 
have  been  successfully  described  w'ith  the  assumption  of  a  constant  flux  layer.  The  state 
variables,  potential  temperature  and  specific  humidity,  obey  the  standard  diffusion  equation 
(Businger,  1973),  perhaps  with  a  source  term  accounting  for  the  moisture  and  temperature  chance 
due  to  injected  or  falling  aerosols  within  an  air  parcel.  Strictly  speaking,  the  conservative 
temperature  variable  is  potential  temperature,  but  within  the  first  hundred  meters  above  the  ocean 
surface,  the  difference  between  potential  temperature  and  the  temperature  is  negligible. 
Henceforth  temperature  will  be  considered  as  the  conservative  heat  variable. 


The  conservation  equation  for  the  state  variable  X  is 


(li 


.  uDi. ,  Codes 
Avail  and/or 
-tst  Speoiai 


□  □ 


MODELLING  THE  F ATE  AND  INFLUENCE  OF  MARINE  SPRa  Y 


where  t  is  time,  z  is  the  vertical  coordinate,  K,  is  the  diffusion  coefficient  for  state  variable  X, 
and  S,  is  the  source  function.  In  this  case  the  source  function  is  the  latent  heat  of  water 
evaporating  or  condensing  for  temperature,  and  the  specific  humidity  for  moisture. 


Aerosol  Moisture  Source  Function 

The  atmospheric  moisture  source  consists  of  the  water  vapor  by  the  evaporating  drop.  The 
amount  of  mass  evaporated  from  the  drop  is  provided  by  Piuppacher  and  Klett  (  1976)  as 

=  4*  CD¥(p.-Pj,(r))  (2) 


where  m  is  condensed  mass  in  gm,  C  is  a  growth  and  ventilation  constant.  Dv  is  the  diffusiviry 
of  water  vapor  (cm:  sec-1),  p_  is  the  ambient  specific  humidity  in  the  region  of  the  particle,  and 
p„,(r)  is  the  saturation  specific  humidity  at  the  surface  of  an  aerosol  of  radius  r. 

The  saturation  vapor  pressure  over  a  droplet  of  size  p,Jr)  is  given  by  the  Koehler  equations 
(Pruppacher  and  Klett,  1978), 


Pj-r)  _  {  A  B 
Pj*)  r  r3 


(3) 
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where  v  is  the  van’t  Hofft  factor,  m,  is  the  mass  of  salt  in  the  droplet,  and  M,  is  the  gram 
molecular  weight  of  salt.  The  Koehler  equation  is  evaluated  with  units  of  radius  in  cm. 
temperature  T  in  Kelvin,  and  mass  in  grams.  The  second  term  is  responsible  for  the  effect  of 
solute  on  evaporation.  In  this  work  we  consider  the  solutes  shown  in  Table  1.  The  solutes  are 


Table  1.  Chemical  constants  of  salts  in  sea  water. 


Salt 

van’t  Hofft  Factor 

Molecular  Weigh: 

NaCl 

2  |  58454 

Na.NO, 

2 

85.005 

(NHj-SO, 

3  j  132.146 

NS 
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considered  to  have  a  uniform  mass  of  10' 15  gm.  This  was  chosen  for  future  consideration  of 
smaller  particles,  since  we  found  the  solute  effects  to  be  negligible  except  for  particles  smaller 
than  10  cm  in  radius.  Testing  is  being  done  on  the  effects  of  different  compositions  on  the 
particle  growth  in  the  context  of  this  modeling  effort. 


Temperature 

The  temperature  change  of  the  system  is  attributed  to  the  flux  gradient  of  temperature  and  a 
source  from  the  latent  heat  absorbed  or  released  by  the  aerosol  particles  within  a  parcel.  The 
general  equation  (1)  becomes 


&  T  _  '  g  8T\  dm 

Tt  "  bz\  HTz)  s~dt 


(5) 


where  the  latent  heat  of  evaporation  of  water  is  given  by  (Pruppacher  and  Klen,  1978)  as 


Ls  =  597.3 


273.— 1  y  =  0.167  ♦  .000367  -7*. 


(6) 


The  Marine  Aerosol  Model 

The  aerosol  population  varies  because  of  local  condensation  and  evaporation,  gravitational 
deposition,  turbulent  transfer,  and  injection  from  the  wave  covered  surface.  The  processes  can 
be  separated  into  continuous  processes,  in  which  mass  is  considered  to  change  continuously,  and 
discrete  processes,  in  which  the  aerosol  number  is  changed  by  injection  or  dynamic  motions.  The 
continuous  processes  are  considered  by  using  the  change  of  a  piecewise  continuous  function; 
discrete  processes  are  described  with  a  conventional  gridded  bin  structure. 

The  particle  population  is  separated  into  a  number  P  of  adjacent  bins.  Within  each  bin  the 
number  density  is  assumed  to  vary  linearly  with  radius.  The  particle  number  distributions  are 
then  given  by, 


(7) 


where  r,  is  the  mid  point  of  the  radius  bin,  and  n,0.k,0  are  determined  from  the  aerosol 
distribution  observable  parameters 

The  number  concentration  of  particles  in  an  air  parcel  with  sices  between  the  limits  of  a  bin  is 
given  by. 
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iV  =  /rj  nrdr 

Jri 

=  ~  Kp  Tl)RX*krJk 


(8) 


where  we  use  the  shorthand  notation, 


The  total  mass  within  a  size  bin  is  gwen  by 


=  Jy  Pr  Wnrfi-krp  r,  )  5  ]• 


Continuous  Transfer  of  Aerosols  -  Evaporation  and  Condensation 
The  change  in  dimension  resulting  from  the  additions  of  mass  is  given  by 

'  3  f  AAfpnpl1'3 
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where  the  assumption  is  made  that  the  change  in  mass  is  distributed  uniformly  over  all  panicles. 
Following  this  change  in  mass,  the  number  of  panicles  transferred  to  another  bin  consists  of  those 
panicles  which  have  grown  (shrunk)  to  a  size  larger  (smaller)  than  the  current  bin  limits.  The 
number  of  panicles  growing  to  the  next  larger  size  is  given  by 

bN  ■=  f'1  n/ir  (12) 


Discrete  Transfer  of  Aerosols  Sedimentation 

The  sedimentation  of  panicles  within  a  layer  was  calculated  assuming  a  uniform  fall  speed  within 
f  ach  size  bin.  Droplet  fall  speeds  were  calculated  using  the  fall  speed  of  a  panicle  of  radius 
equal  to  the  mean  panicle  radius  in  the  bin.  The  fall  speed  was  the  unmodified  Stokes  fall  speed, 


2  ga 
9v 


(13) 


The  change  in  the  number  of  panicles  in  the  bin  is  given  by 


AN"  =  N,"F,  -  Aft  F 


(14) 


where  N  are  the  panicle  number  concentrations,  and  F  are  fractions  of  panicles  lost  or  gained 
by  sedimentation.  The  fraction  of  panicles  leaving  a  layer  during  a  time  interval  is 


F 


vA T 
A  z 


(15) 


The  time  interval  and  mean  size  must  be  chosen  judiciously  so  that  panicles  do  not  uniformly 
fall  through  their  air  parcels. 

Injection 

The  aerosol  source  function  is  considered  to  be  an  actively  breaking  wave,  providing  an  aerosol 
source  injecting  panicles  with  a  range  of  upward  venical  velocities.  Rather  than  an  Eulerian 
scheme  tracking  each  panicle,  we  model  the  phenomena  with  an  injection  functionwhich  is  height 
dependent.  The  source  function  at  the  surface  is  the  source  function  of  Marks  and  Nlcnahan 
(1988).  This  function  is  translated  to  the  parameters  of  a  modified  Gamma,  or 
Nukiyama-Tanasawa  distribution,  using  a  maximum  likelihood  method  for  choice  of  parameter^ 
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(Goroch,  1980).  This  function  is  dependent  on  wind  speed  with  a  power  law  dependence, 
following  observations  of  Monahan  et  al.  (1983).  The  vertical  distribution  is  exponential,  to 
reproduce  observations  of  DeLeeuw  (1986). 

Initial  Conditions 

The  initial  conditions  of  the  model  were  calculated  assuming  a  well  mixed  atmospheric  surface 
layer,  defined  by  the  10m  value  of  temperature,  specific  humidity  and  wind  speed.  The  vertical 
profiles  of  these  were  initially  calculated  using  Monin-Obukhov  scaling  as  reviewed  by  Smith 
(1989).  The  three  relevant  equations  for  wind  speed,  specific  humidity  and  potential  temperature 
are,  respectively: 


(16) 


(17) 


The  flux  quantities  were  obtained  from  the  initial  conditions  by  using  the  flux  profile  calculations 
of  Smith  (1988). 

Summary 

This  model  of  the  aerosol  dynamics  in  the  marine  boundary  layer  emphasizes  the  microphysics 
of  particles  in  the  marine  boundary  layer.  The  several  processes  affecting  the  motion  and  growth 
of  particles  are  separated  into  larger  scale  changes  of  the  temperature  and  moisture,  and  the 
discrete  and  continuous  change  of  the  aerosol  number  and  mass  distributions.  The  method  is 
applied  to  a  Lagrangian  model,  but  could  as  well  be  applied  to  an  Eulerian  description  of  aerc  :>1 
dynamics. 


AX.  GOROCH  A  CONTINUOUS  BIN  MODEL  OF  AEROSOL  BEHAVIOR  IN  THEM  ARISE  BOUNDARY  LAYER 


References 


Bortkovskii,  R.S.,  1987:  Air-Sea  Exchange  of  Heat  and  Moisture  During  Storms.  Revised 
English  Edition,  Edward  C.  Monahan,  ed.,  D. Seidel  Publishing,  Dordrecht. 

Burk,  Stephen  D.,  1984:  "The  Generation,  Turbulent  Transfer  and  Deposition  of  the  Sea-Salt 
Aerosol,"  J.  Atmos.  Sci..  41  (20),  3040-3051. 

Businger,  Joost,  1973:  "Turbulent  Transfer  in  the  Atmospheric  Surface  Layer,"  in  Workshop  on 
Micrometeorology.  Amer.  Meteor.  Soc.,  Boston,  MA. 

De  Leeuw,  G„  1986:  "Vertical  profiles  of  giant  panicles  close  above  the  sea  surface".  Telius. 
38B.  51-61. 

Fairall,  C.  W.,  and  J.  B.  Edson,  1989:  "Modeling  the  Droplet  Contribution  to  the  sea-to-air 
Moisture  Flux,"  in  Climate  and  Health  Implications  of  Bubble-mediated  Sea-Air  Exchange, 
Edward  C.  Monahan  and  Margaret  Van  Patten,  eds.,  distributed  by  the  Connecticut  Sea  Grant 
College  Program,  Marine  Sciences  Institute,  University  of  Connecticut  at  Avery  Point,  Groton, 
CT  06340. 

Goroch,  Andreas  K.,  1980:  "Optimal  Fitting  of  Marine  Aerosol  Distributions  Proc.  2nd  Conf. 
Coastal  Meteorol.,  Los  Angeles,  CA  p  136-140. 

Ling,  S.C.,  and  T.W.  Kao,  1976:  "Parameterization  of  the  moisture  and  heat  transfer  process  over 
the  ocean  under  whitecap  sea  states".  J.  Phvs.  Oceanogr..  6,  306-315. 

Marks,  Roland  and  Edward  C.  Monahan,  1988:  "Relationships  Between  Marine  Aerosols, 
Oceanic  whitecaps,  and  Low-Elevation  winds  Observed  During  the  HEXMAX  Experiment  in  the 
North  Sea,"  in  Climate  and  Health  Implications  of  Bubble-mediated  Sea-Air  Exchange.  Edward 
C.  Monahan  and  Margaret  Van  Patten,  eds.,  distributed  by  the  Connecticut  Sea  Grant  College 
Program,  Marine  Sciences  Institute,  University  of  Connecticut  at  Avery  Point,  Groton,  CT  06340. 

Monahan.  E.C.,  D.  E.  Spiel  and  K.L.  Davjdson.  1986:  "A  Model  of  Marine  Aerosol  Generation 
via  Whitecaps  and  Wave  Disruption,"  in  Oeanic  Whitecaps.  D.  Reidel.  Dordrecht,  Holland,  pp 
167-174. 

Pruppacher,  H.  R.  and  J.O.  Klett,  1978:  Microphysics  of  Clouds  and  Precipitation.  Reidel, 
Dordrecht,  Holland. 

Smith,  Stuart  D.,  1988.  "Coefficients  for  Sea  Surface  wind  Stress.  Heat  Flux,  and  Wind 
Profiles  as  a  Function  of  Wind  Speed  and  Temperature,"  J,  Geophvs.  Res.  93  (02), 
15467-15472. 


153 


* 


MODELLING  THE  FATE  AND  INFLUENCE  OF  MARINE  SPRAY 


Stramska,  M.,  1986:"  Vertical  Profiles  of  sea-salt  aerosol  in  the  atmosphere  surface  layer  A 
Numerical  Model".  Preprint,  to  appear  in  Acta  Geophvs.  Pol. 

Young,  Kenneth  C.,  1974:  "A  Numerical  Simulation  of  wintt.nime  Orographic  Precipitation:  Part 
I.  Description  of  Model  Microphvsics  and  Numerical  Techniques."  J.  Atmos.  Sci.  3 1 . 1735-1748. 


NOARL  contribution  no.  90:075:442.  Approved  for  public  release;  distribution  is  unlimited. 


154 


